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We review the developments of QCD multipole expansion and its applications to hadronic tran- 
sitions and some radiative decays of heavy quarkonia. Theoretical predictions are compared with 
updated experimental results. 

I. INTRODUCTION 



Heavy quarkonia are the simplest objects for studying the physics of hadrons due to their nonrelativistic nature. 
Although the spectra of heavy quarkonium systems cc and bb have been successfully explained by certain QCD 
motivated potential models, some of their decays concerning nonperturbative QCD are difficult to deal with. 
Hadronic transitions 
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y—{ ' are of this kind. In (1), <£>/, $>p and h stand for the initial state quarkonium, the final state quarkonium, and the 
emitted light hadron(s), respectively. Hadronic transitions are important decay modes of heavy quarkonia. For 
instance, the branching ratio for ip 1 — > J/tp + n + n is approximately 50%. 

In the cc and bb systems, the typical mass difference — M$ p is around a few hundred MeV, so that the 
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typical momentum of the light hadron(s) h is low. So far as the coupled-channel effect is not concerned, the 
light hadron(s) h are converted from the gluons emitted by the heavy quark Q or antiquark Q in the transition. 
[ So the typical momentum of the emitted gluons is also low, and thus perturbative QCD does not work in these 
' processes. Certain nonperturbative approaches are thus needed for studying hadronic transitions. In this article, 
we review the theoretical framework and applications of a feasible approach, QCD multipole expansion (QCDME), 
which works quite well in predicting hadronic transition rates in the cc and bb systems. In addition to hadronic 
transitions, QCDME can also lead to successful results in certain radiative decay processes such as J/ip — ► 77; 
and J/tp — > 777'. 

This paper is organized as follows. In Sec. II, we review the theoretical framework and the formulation of 
QCDME. Sec. Ill deals with applications of QCDME to various hadronic transition processes in the nonrelativistic 
single-channel approach including hadronic transitions between S'-wave quarkonia, between P-wave quarkonia, 7T7t 
transition of the D-w&ve quarkonia, and the search for the spin-singlet P-wave quarkonium h c through hadronic 
transition. Sec. IV is on the nonrelativistic coupled-channel theory of hadronic transitions. In Sec, V, we show 
how QCDME makes successful predictions for the radiative decays J/ip — > 777 and J/ijj — > 7?/, etc. A summary 
is given in Sec. VI. 



II. QCD MULTIPOLE EXPANSION 

Multipole expansion in electrodynamics has been widely used for studying radiation processes in which the 
electromagnetic field is radiated from local sources. If the radius a of a local source is smaller than the wave 
length A of the radiated electromagnetic field such that a/ A ~ ak < 1 (k stands for the momentum of the photon), 
ak can be a good expansion parameter, i.e., we can expand the electromagnetic field in powers of ak. This is the 
well-known multipole expansion. In classical electrodynamics, the coefficient of the (ak) term in the multipole 
expansion contains a factor — . Hence multipole expansion actually works better than what is expected 
by simply estimating the size of (ak) \ 



Due to the nonrelativistic nature of heavy quarkonia, the bound states of a heavy quark Q and its antiquark 
Q can be calculated by solving the Schrodinger equation with a given potential model, and the bound states are 
labelled by the principal quantum number n, the orbital angular momentum L, the total angular momentum J, 
and the spin multiplicity a (u = 1 or 3), i.e., n a Lj. The typical radius a = \J (r 2 ) of the cc and bb quarkonia 
obtained in this way is of the order of 1CT 1 fm. With such a small radius, the idea of multipole radiation can be 
applied to the soft gluon emissions in hadronic transitions. Consider an emitted gluon with a momentum k. For 
typical hadronic transition processes, k ~ few hundred MeV, so that ak is of the order of 10 _1 . Thus multipole 
expansion works for hadronic transition processes. Note that the convergence of QCDME does not depend on 
the value of the QCD coupling constant g s . Therefore QCDME is a feasible approach to the soft gluon emissions 
in hadronic transitions (1). 

QCDME has been studied by many authors [1-6]. The gauge invariant formulation is given in Ref. [5]. Let 
ip(x) and Ap(x) be the quark and gluon fields, respectively. Following Refs. [5], we introduce 



^A;'(x,t) = U-\x,t)^-A a Jx)U{x,t) - -U- 1 (x,t)d f ,U(x,t), 
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where U(x,t) is defined by [6] 



U(x,t) = Pexp 



^A a (x',t)-dx' 



(2) 



(3) 



in which P is the path-ordering operation, the line integral is along the straight-line segment connecting the 
two ends, X = (x\ + x 2 )/2 is the center of mass position of Q and Q, and x denotes x± or x 2 . With these 
transformed fields, the part of the QCD Lagrangian related to the heavy quarks becomes [5] 
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^{x 2 ,t) 1 ^{x 2l t)d 3 x 1 d 3 x 2 , 
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where \q/2 = 1. Note that the transformed quark field ^(a;, t) is dressed with gluons through U^ 1 (x^ t) defined 
in (3). We see from Eq. (4) that the dressed quark field ^(a^t) serves as the constituent quark field interacting 
via the static Coulomb potential in the potential model. In addition, it is the transformed gluon field A^ (not 
the original A°j) that appears in the covariant derivative in (4). A a ' contains non-Abelian contributions through 
U(x,t). 

Following Ref. [5] , we generalize the Coulomb potential in Eq. (4) to the static potential including the confining 
potential in potential models, and we write down the following effective Lagrangian [5] 
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J a=0 
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^(x 2 ,t) 1 a ^(x 2 ,t)d 3 x 1 d 3 x 2 , 
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where Vi(|aJi — x 2 \) is the static potential (including the confining potential) between Q and Q in the color-singlct 
state, and V^(|a3i — x 2 \) is the static potential between Q and Q in the color-octet state. This Cq relates the 
QCD Lagrangian to the potential models. 

Now we consider the multipole expansion. Inside the quarkonium, \x — X\ < a. So we can make an expansion 
by expanding the gluon field A^(x,t) in Taylor scries of x — X at the center of mass position X. The Taylor 
series is an expansion in powers of the operators (x — X) • V and (x — X) X V applying to the gluon field. 
After operating on the gluon field with the gluon momentum k, these operators are of the order of ak. This is 
QCDME. It has been shown in Ref. [5] that this operation leads to 



A$(x, t) = A a '{X, t)-{x-X)- E a (X, t) + ■ 
A a '(x,t) = -l(x-X) x B a (X,t) + ---, 



(6) 
(7) 



where E a and B a are color-electric and color-magnetic fields, respectively. 

In Ref. [5], the corresponding Hamiltonian was derived based on the above formulation. This is more convenient 
in using the nonrelativistic perturbation theory. The obtain Hamiltonian is [5] 

rreff _ rr(0) tt(1) 

n QCD ~ n QCD ' n QCDi \°) 

where 



H qcd = j ^(x u t)^{x u t)H^(x 2 ,t)^(x 2 ,t)d 3 Xid 3 x 2 



(9) 



with 



and 



H = -^(df + d 2 2 + Vi(|xi - x 2 \) + £ y y^ 2 (|xi - x 2 |) + 2m Q , (10) 

^ a— 1 



Hqcd = Hi+H 2l 

Hi=Q a A%(X,t), H 2 = -d a -E a (X,t)-m a -B a (X,t) + ---, (11) 

in which 

Q a = g E j * t) y*( x > ( 12 ) 

d a = g E f(x-X)^(x,t)^(x,t)d 3 x, (13) 

m a = ^J(x-X)x *t( X) t) 7 ^*(x, i)<2 3 x (14) 

are the color charge, color-electric dipole moment, and color-magnetic dipole moment of the QQ system, respec- 
tively. Note that Eq. (5) is regarded as an effective Lagrangian. Considering that the heavy quark may have an 
anomalous magnetic moment, we have taken in Eqs. (12), (13) and (14) the symbols g E and qm to denote the 
effective coupling constants for the electric and magnetic multipole gluon emissions, respectively. We shall see 
later in Sec. Ill that taking a E and au as two parameters is needed phcnomenologically. 

We are going to take Hq CD as the zeroth order Hamiltonian, and take Hqcd as a perturbation. This is 
different from the ordinary perturbation theory since Hq^ cd is not a free field Hamiltonian. Hq CD contains 
strong interactions in the potentials in H, so that the eigenstates of Hq CD arc bound states rather than free 
field states. For a given potential model, the zeroth order solution can be obtained by solving the Schrodinger 
equation with the given potential. Moreover, we see from Eqs. (12), (13) and (14) that only H 2 in H Q 1] CD is of 
0(ak), while Hi is of 0((ak) ). So that we should keep all orders of Hi in the perturbation expansion. 

The general formula for the S matrix element between the initial state \I) and the final state \F) in this 
expansion has been given in Ref. [6] , which is 



(F\S\I) = -i2n5(E F +uo F - E t )(f 



H 2 j- H 2 ■ ■ ■ H 2 -t H 2 

E i - h qcd + ido -Hi Ej — H Q % D + id» - Hi 
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where lo e is the energy of the emitted gluons. This is the basis of the study of hadronic transitions in QCDME. 
Explicit evaluation of the S matrix elements in various cases will be presented in Sec. III. 



III. PREDICTIONS FOR HADRONIC TRANSITIONS IN THE SINGLE-CHANNEL APPROACH 

In this section, we shall show the predictions for hadronic transitions rates in the single-channel approach 
(inclusion of coupled-channel contributions will be given in Sec. IV) . In this approach, the amplitude of hadronic 
transitions (1) is diagrammatically shown in FIG. 1 in which there are two complicated vertices: namely, the vertex 
of multipole gluon emissions (MGE) from the heavy quarks and the vertex of hadronization (H) describing the 
conversion of the emitted gluons into light hadron(s). The MGE vertex is at the scale of the heavy quarkonium, 
and it depends on the property of the heavy quarkonium. The H vertex is at the scale of the light hadron(s), 
and is independent of the property of the heavy quarkonium. In the following, we shall treat them separately. 



h 




FIG. 1: Diagram for a typical hadronic transition in the single-channel approach. 



Hadronic Transitions Between 5- Wave Quarkonia 



Let us first consider the case of tttt transitions between S'-wave quarkonia, n^Si — > npSi + it + n. These 
processes are dominated by double electric-dipole transitions (E1E1). The transition amplitude can be obtained 
from the S matrix element (15). With certain algebra, we obtain [5-7] 
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where x is the separation between Q and Q, and (D )bc = dbedo — g s fabcA^. Let us insert a complete set of 
intermediate states with the principal quantum number K and the orbital angular momentum L. Then Eq. (16) 
can be written as 
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(17) 



According to the angular momentum selection rule, the intermediate states must have L = L' = 1. The 
intermediate states in the hadronic transition are the states after the emission of the first gluon and before the 
emission of the second gluon shown in FIG. 1, i.e. they are states with a gluon and a color-octet QQ. There is 
strong interaction between the gluon and the color-octet QQ since they all carry colors. Thus these states are 
the so-called hybrid states. It is difficult to calculate these hybrid states from the first principles of QCD. So 
we shall take a reasonable model for it. The model should reasonably reflect the main properties of the hybrid 
states and should contain as few unknown parameters as possible in order not to affect the predictive power of the 
theory. There is a quark confining string (QCS) model [8] satisfying these requirements. The QCS model is a one 
dimensional string model in which the strong confining force between Q and Q is described by the ground-state 
string, and gluon excitation effects are described by the vibrations of the string [8] . Our intermediate states QQg 
are thus described by the first vibrational mode in this model. The QCS model is not the only one satisfying the 
above requirements. Another possible model satisfying the requirements is the MIT bag model for the hybrid 
states. Hadronic transitions with the MIT bag model as the model for the intermediate states has been studied 
in Ref. [9]. It is shown that, with the same input data, the predictions in this model are very close to those in the 
QCS model although the absolute intermediate state energy eigenvalues in the bag model are much higher than 
those in the QCS model. Thus the predictions are not sensitive to the specific energy spectrum of the intermediate 
states. In the following, we take the calculations with the QCS model as examples. Explicit calculations of the 
first vibrational mode as the intermediate states is given in Ref. [7]. With this model, the transition amplitude 
(17) becomes [7] 



g\ ^ ^ F \x k \KL){KL\x^i) ■ 
Meiei = l -^2-^ El _ E[(L \ n7T \ E k E l 1°). 



(18) 



where Ekl is the energy eigenvalue of the intermediate vibrational state \KL). We see that, in this approach, 
the transition amplitude contains two factors: namely, the heavy quark MGE factor (the summation) and the H 
factor (7T7r|i^_E ; a |0). The first factor concerns the wave functions and energy eigenvalues of the initial and final 



state quarkonia and the intermediate states. These can be calculated for a given potential model. Let us now 
consider the treatment of the second factor. The scale of the H factor is the scale of light hadrons which is very 
low. Therefore the calculation of this matrix element is highly nonperturbative. So far, there is no reliable way 
of calculating this H factor from the first principles of QCD. Therefore we take a phenomenological approach 
based on an analysis of the structure of this matrix element using PCAC and soft pion technique in Ref. [10]. 
In the center-of-mass frame, the two pion momenta qi and q 2 are the only independent variables describing this 
matrix element. According Ref. [10], we can write this matrix element as [7] 



?f(n a ( qi )n ( q2 )\E a k E?\O) = J aP 

V(2wi)(2w 2 ) 



CiSkiQi + C 2 ^<zifc<72; + qimk - ^Suqi ■ q^j 



(19) 



where C\ and C 2 are two unknown constants. In the rest frame of <&/ (the center-of-mass frame), for a given tttt 
invariant mass M WTr , the C\ term is isotropic (5- wave), while the C 2 term is angular dependent (D-wave). In 
the nonrelativistic single-channel approach, the MGE factor in (18) is proportional to Ski due to orbital angular 
momentum conservation. So that only the C\ term contributes to the 5-state to 5-state transitions. In this case, 
the n|5i — > n F S\ + n + n transition rate can be expressed as [7] 

r(nfS! - n 3 ^ tt tt) = |Ci| 2 G|0, F ol 2 . (20) 
where the phase-space factor G is [7] 
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in which Rj, Rp, and Rkl are radial wave functions of the initial, final, and intermediate vibrational states, 
respectively. These radial wave functions are calculated from the Schrodinger equation with a given potential 
model. 

Now there is only one overall unknown constant C\ left in this transition amplitude, and it can be determined 
by taking a well-measured hadronic transition rate as an input. So far, the best measured 5-state to S'-state tttt 
transition rate is Y{ip' — > J/ip nn). The updated experimental value is [11] 

r tot (V/) =281±17keV, 

B(lp' -» J/l/> TT+TT-) = (31.8 ± 1.1)%, 

B(ip' -» J /il) 7r°7r°) = (18.8 ± 1.2)%. (24) 

We take this as an input to determine C\. Then we can predict all the 5-state to 5-state tttt transitions rates 
in the T system. Since the transition rate (20) depends on the potential model through the amplitude (23), 
the determined value of \C\\ is model dependent. In the following, we take the Cornell Coulomb plus linear 
potential model [12] and the Buchmuller-Grunberg-Tye (BGT) potential model [13] as examples to show the 
determined |Ci| and the predicted rates of T' — > T tttt, T" — ► T tttt, and T" — > T' tttt. The results are listed 
in TABLE I [59]. We see that the predicted ratios T(T" -> T 7r7r)/L(T' — > T tttt) w 1.2/7.8 = 0.15 and 
r(T" — » T' tttt)/T(T' — > T tttt) w 0.53/7.8 = 0.07 in the BGT model are close to the corresponding experimental 
values L(T" -» T 7r7r)/L(T' -» T tttt) w 1.72/12.0 = 0.14 and L(T" -» T' 7T7r)/L(T' -» T tttt) w 1.26/12.0 = 
0.11. However, the predicted absolute partial widths are smaller than the corresponding experimental values by 
roughly a factor of (50— 75)%. Moreover, when the distributions are considered, the situation will be more 
complicated. We shall deal with these issues in Sec. IV. 



TABLE I: The determined |C*i| 2 and the predicted rates T(T' ->■ T tttt), T(T" ->■ T tttt), and T(T" -> T' tttt) ( in kcV) in the 
Cornell model and the BGT model. The corresponding updated experimental values of the transition rates quoted from Ref. [11] 
are also listed for comparison. 



Cornell BGT Expt. 



|Ci| 2 83.4 x 1(T 6 67.8 x 1(T 6 

T(T' -f Ttttt) (keV) 8.6 7.8 12.0 ± 1.8 

r(T" -f Ttttt) (kcV) 0.44 1.2 1.72 ± 0.35 

T(T" -> T'tttt) (keV) 0.78 0.53 1.26 ± 0.40 



Note that the phase space factor G in T" — > T 7T7t is much larger than that in T' — > T 7T7t, G(T" — > 
T 7T7r)/G(T' — » T 7T7r) = 33. So one may naively expect that T(T" — > T 7T7t) > T(T' — > T 7T7t). However, 
as we see from the experimental values in TABLE I that T(T" -> T 7T7r)/r(T' -> T tttt) w 1.72/12.2 = 0.14. 
The reason why our predictions for this ratio is close to the experimental value is that the contributions from 
various intermediate states to the overlapping integrals in the summation in f^l Q [cf. Eq. (23)] drastically cancel 
each other due to the fact that the T" wave function contains two nodes. This is a characteristic of this kind of 
intermediate state models (QCS or bag model). To see this, let us take a simplified model for the intermediate 
states and look at its prediction. If we make a simplification assumption that the variation of the factor 1/(2?/ — 
Hqq D — iDo) in Eq. (17) is sufficiently slow such that it can be approximately represented by a constant which 
can be taken out of the summation, the summations and in (17) can then be carried out and the 

KL K'L> 

double overlapping integrals in the numerator in (23) reduces to a single integration j R^,, (r)r 2 2?x (r)r 2 dr . This 
simplified model predicts a rate T(T" — > T tttt) larger than the experimental value by orders of magnitude [7]. 
Therefore, taking a reasonable model for the intermediate states is crucial for obtaining successful predictions in 
the QCDME approach to the MGE factor. 

The transitions nf Si — > ripSi + r\ are contributed by E1M2 and M1M1 transitions, and is dominated by the 
E1M2 transition. The transition amplitude is 

^^=-2^—2^ Ej-E KL H2? fc 9 m ^|0), (25) 

where S is the total spin of the quarkonium. The MGE matrix element is proportional to 5km- Similar to the 
idea in (19), we can phenomenologically parameterize the hadronization factor according to its Lorentz structure 
as [14] 

i^( v ( q )\E a k d k Bf\0) = z(2tt) 3 / 2 C 3 qi , (26) 
6 

in which the phenomenological constant can be determined by taking the data [11] 

r tot (^') = 277 ±22 keV, B(ip' -» J/ip rj) = (3.17 ± 0.21)% (27) 
as input, and so we can predict the rates for T' — > T r\ and T" — > T 77. This is equivalent to 



r(T(n?Si) - T V ) = 





2 




\q(bb)\ 3 



/201o( cc ) 



l<?(cc)| ; 



Tty' J/V, V ). (28) 



where q(bb) and q(cc) are the momenta of T] in T(n^Si) — > T i] and -0' — * J ftp V: respectively. Taking the BGT 
model as an example to calculate the ratio of transition amplitudes in (28), we obtained 

T(T' -> T ry) = 0.022 keV, T(T" -> T r?) = 0.011 keV. (29) 

These are consistent with the present experimental bounds [11] 

T(T' -» T 77) < 0.086 kcV, T(T" -» T' 77) < 0.058 kcV. (30) 



We can also compare the ratios R' = T(T' — > T f])/T(^' -» J/V> r?) and P" = T(T" -» T rf)/T{^' -> J /if) rf) with 
the recent experimental measurements. Recently BES has obtained an accurate measurement of T(ip' — > J/^ 77) 
and r(V>' -> J/> tt°) [15]. With the new BES data and the bounds on T(T' -> T 77) and T(T" -» T 77) [11], the 
experimental bounds on R' and J?" are [15] 

R'Uxpt < 0.0098, R"\ expt < 0.0065. (31) 

Taking the BGT model to calculate the ratios R' and R", we obtain 

R'\bgt — 0.0025, R"\bgt = 0.0013. (32) 

These are consistent with the new experimental bounds (31). 



B. 7777 Transitions Between P-Wave Quarkonia 



Let us consider the hadronic transitions 2 3 P/ J — > l 3 Pj F + ir + ir. For simplicity, we use the symbol T(Jj — > Jp) 
to denote r(2 3 Pj, -» l 3 Pj F 7r7r). These are also dominated by E1E1 transitions. The obtained results are [5, 7] 



r(o-o) = l|c 1 | 2 G 



r(o -» 2) = 5r(2 -» 0) 
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r(o -» 1) = r(i -» 0) = 0, 
2 
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r(i - i) = r(o - o) + -r(o ^ 2), 



r(i ^ 2) = ^r(2 1) = ^r(o 2), r(2 ^ 2) = r(o - 0) + ^r(o 2), 



(33) 



where the phase-space factor H is 
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20 Ms, 



2 /f 2 
1 + ^ 



3M. 



(M^ + 2m 2 M 2 T + 6m4) 



dM 2 T ,(34) 



with X defined in Eq. (22). 

Now the rates in (33) depend on both C\ and C^. We know that C\ has been determined by the input (24). So 
far, there is no well measured hadronic transition rate available for determining the ratio C2/C1. At present, to 
make predictions, we can only take certain approximation to estimate CijC\ theoretically. The approximation 
taken in Rcf. [7] is to assume that the H factor (irir\E%E®\0) can be approximately expressed as 



(™\E%E?\0)cc{gg\E%E?\0), 



(35) 



i.e., (7T7r|-E£-E"|0) approximately contains a factor {gg\E%Ef\0) and another factor describing the conversion of 
the two gluons into 7T7t which is assumed to be approximately independent of the pion momenta in the hadronic 
transitions under consideration. The R.H.S. of Eq. (35) can be easily calculated. Comparing the obtained result 
with the form (19), we obtain 



C2/C1 



(36) 



in such an approximation. This is a crude approximation which can only be regarded as an order of magnitude 
estimate. So it is likely that C2/O1 <~ 0(1) rather than O(10 _1 ) or O(10). A reasonable range of O2/O1 is 



1 < O2/O1 < 3. 



(37) 



With this range of O2/O1, the obtained transition rates T(Jj — > Jp) of Xf ) (2 3 Pj J ) — ► Xb{l 3 P.j F )irir in the 
Cornell model [12] and the BGTmodel [13] are listed in TABLE II. The relations between different T(J/ -> J F ) 
given in Eq. (33) reflect the symmetry in the E1E1 multipole expansion [5], so that experimental tests of these 
relations are of special interest. Recently, CLEO reported a preliminary observation of the hadronic transitions 
Xb(2 3 Pj 7 ) — > Xb(l 3 P.j F ) 7T7T f° r Ji = Jf = 1 an d 2 [16]. In a very recent paper [17], CLEO measured the transition 
rate, and the obtained result is r(x 6 (2 3 P 7l ) -> Xb(^Pj F )^) = (0.83 ± 0.22 ± 0.08 ± 0.19) keV (J/ = J F = 1, 2) 
which is consistent with the predicted rates T(l — > 1), and T(2 — » 2) listed in TABLE II. 



TABLE II: Predicted transition rates F(Ji — > Jf) of Xf>(2 3 -Pj 7 ) — > Xf>(l 3 -P./F) 7r7r with the parameter range (37) in the 
Cornell model [12] and the BGTmodel [13]. 



Model 






r(Jj -> Jf) (keV) 










r(o -> o) 


r(o -> 2) 


r(i-i) 


r(i- 


-2) 


r(2 -»■ 2) 


Cornell 


0.4 


0.004-0.04 


0.4 


0.003 


-0.03 


0.4 


BGT 


0.4 


0.002-0.02 


0.4 


0.001 


-0.01 


0.4 



C. 7T7T Transitions of D-Wave Quarkonia 

■0(3770) (or 0") is commonly regarded as essentially the ID state of the charmonium. It lies above the DD 
threshold, so that it is usually believed that 0(3770) mainly decays into the open channel DD. Experimental 
observations show that the directly measured -0(3770) production cross section at e+e~ colliders is [18, 19] 

ct(0(377O)) = 7.5 ±0.8 nb, (38) 

while the e+e" — > 0(3770) — > DD cross section is [20] 

ct(0(377O) -> DD) = 5.0 ± 0.5 nb. (39) 

This discrepancy may indicate that there are considerable non-DD decay modes of 0(3770). One of the possible 
non-DD decay modes is the hadronic transition 0(3770) — > J/ip tttt. Theoretical studies of hadronic transitions of 
the _D-wave quarkonia have been carried out by several authors in different approaches leading to quite different 
predictions [7, 21-24]. In the following, we briefly review the approach given in Rcfs. [23, 24], and compare the 
predictions with the recent experimental result and with other approaches. 

The measured leptonic width of 0(3770) is (0.26 ± 0.04) keV [11]. If we simply regard 0(3770) as a pure ID 
state of charmonium, the predicted leptonic width will be smaller than the experimental value by an order of 
magnitude. Therefore people consider 0(3770) as a mixture of charmonium states [23-25]. State mixing is an 
important consequence of the coupled-channel theory, especially for states close to or beyond the open channel 
threshold. Take a successful coupled-channel model, the unitary quark model (UQM) [26], as an example. In 
this model, 0(3770) is a mixture of many 5- wave and D-wave states of charmonium, but the main ingredients 
are the tp(lD) and 0(25) states. Neglecting the small ingredients, we can write 0' and 0(3770) as 

0' = 0(25) cos + 0(LD) sin0, 

0(3770) = -0(25) sin 9 + 0(1L>) cos 9. (40) 

The UQM gives 9 ~ — 8° [26]. Instead of taking a specific coupled-channel model, we take a phenomenological 
approach determining the mixing angle 9 by fitting the ratio of the leptonic width of 0' and 0(3770). The leptonic 
widths of 0(25) and 0(1D) are proportional to the wave function at the origin 02s(O) and the second derivative 

of the wave function at the origin -^= - ^' 1 - p (^/^ r respectively. Therefore the determination of 9 depends on 

V2 2mi 

the potential model. Here we take two potential models as illustration: namely, the Cornell potential model [12] 
and the improved QCD motivated potential model by Chen and Kuang (CK) [27] which leads to more successful 
phenomenological results. The determined values of 9 are 

Cornell : = -10°, 

CK : 9 =-12°. (41) 

These are all consistent with the UQM value. There can also be an alternative solution with 9 ~ 30°, but it is 
ruled out by the measured M W7r distribution of 0' — > J /if) + tt + tt. 

This transition is also dominated by E1E1 gluon emission. The transition rate is [23] 



T(0(377O) -> J/0 tttt) = |d | 



sin 2 tfG(0') |.f 2 1 o 1 1 1 o(0')| 2 + ^ 



Ci 



2 



cos 2 iJ(0") |/ 1 1 2 1 1 1 o (0")| 2 



(42) 



This transition rate depends on the potential model through the amplitudes /201b, /mo an d the value of 
C2/C1. We take the Cornell model [12] and the CK model [27] as examples. Taking the possible range for C2/C1 
given in (37), we obtain the values of F(V>(3770) -> J/ip + ir+ + tt") listed listed in TABLE III [60]. Note that 

TABLE III: The predicted transition rate r(i/>(3770) -> J/ip + vr+ + tt~ ) (in kcV) in the Cornell model and the CK model with the 
updated input data (24). 



Model 


L(V>(3770) -f J/ip tt+tt-) (keV) 


Cornell 


26 - 139 


CK 


32 - 147 



S-D mixing only affects a few percent of the rate, so that the rate is essentially T(ip(lD) — > J/ip tt + tt ). 

Recently, BES has measured the rate L(>(3770) -> J/ip + tt+ + tt") based on 27.7 pb" 1 data of ^(3770). The 
measured branching ratio is [28] 

5(^(3770) -> J/ip + 7T+ + tt") = (0.34 ± 0.14 ± 0.09)%. (43) 

With the total width [11] 

L tot 0(3770)) = 23.6 ± 2.7 MeV, (44) 

the partial width is [28] 

T BES (ip{3770) -> J/V> + 7T+ + tt") = 80 ± 32 ± 21 keV. (45) 

This is in agreement with the theoretical predictions in TABLE III. Taking the BES data (45) and Eq. (42) to 
determine C2/C1 , we obtain 

C2/C1 = 2t°-j. (46) 

This shows that C2/C1 is really of O(l). 

Very recently, CLEO-c also detected the channel ?/>(3770) — > J/f/i + 7r + + 7r~ with higher precision, and the 
measured branching ratio is [29] 

#0(3770) -» J/V> + 7T+ + tt") = (0.214 ± 0.025 ± 0.022)%. (47) 

With the V(3770) total width (44), the partial width is 

r(V>(3770) -» J/r/" + 7T+ + tt") = 50.5 ± 16.9 kcV. (48) 

We can also determine C2/C1 from (48) and (42), and the result is 

C2/C1 = 1.52t8;» (49) 

This is consistent with the value (46) determined from the BES data, but with higher precision. 

An alternative way of calculating this kind of transition rate taking the approach to the H factor proposed 
by Ref. [4] was carried out in Ref. [22]. The so obtained transition rate is smaller than the above theoretical 
prediction by two orders of magnitude. So it strongly disagrees with (45) and (48). Therefore the approach given 
in Ref. [4] is ruled out by the BES and CLEO-c experiments. 

For the T system, state mixings are much smaller [26]. Neglecting state mixings, the T(1D) — > T(lS')+7r + +7r~ 
transition rate is proportional to [C2/C1] 2 . Taking the determined values of C2/C1 in (46) and (49), we obtain 
the corresponding transition rates: 1.3 kcV < F(T(1Z?) — > T(15 f ) + ir + + 7r~) < 14 keV [from (46)] and 
2.0 keV < r(T(LD) -» T(1S) + tt+ + tt") < 5.0 keV [from (49)], respectively. The lower values in these ranges 
are consistent with the CLEO bound [30]. Improved measurement of the T(1D) — > T(15) + 7r+ + n~ rate is 
desired. 



D. Searching for the h c States 



The spin-singlet P-wave states (I 1 Pi) are of special interest since the difference between the mass of the 
l l P\ state and the center-of-gravity of the l 3 Pj states M c . og = (5M]3p 2 + 3Misp 1 + Mi3p )/9 gives useful 
information about the spin-dependent interactions between the heavy quark and antiquark. There have been 
various experiments searching for the h c { , 4>(l 1 Pi)) state. 

In the pp collision, h c can be directly produced. In 1992, the E760 Collaboration claimed seeing a significant 
enhancement in pp — > J/ip + 7r° at yfs = 3526.2 MeV which was supposed to be a candidate of h c [31]. However, 
such an enhancement has not been confirmed by the successive E835 experiment from a careful scan in this region 
with significantly higher statistics [32]. Instead, the E835 experiment recently found the h c state via another 
channel pp — > h c — > t? c 7, and the measured resonance mass is M^ c = 3525.8 ±0.2 ±0.2 MeV with a width < 1 
MeV [32]. The measured production rate is consistent with the theoretical range given in Ref. [14] (see Ref. [32]). 

At the e + e~ colliders, the h c state cannot be produced directly in the s-channcl due to its CP quantum 
number. Because of the limited phase space, the best way of searching for the h c state at CLEO-c or BES is 
through the isospin violating hadronic transition [14, 24, 33] 

V>' -^h c + ir°. (50) 

Theoretical calculations of this transition rate considering S-D mixing in ip' and suggestions for tagging the h c 
are given in Ref. [24]. Here we give a brief review of it. 

The process ip' — > h c + 7r° is dominated by E1M1 transition. The transition amplitude is 

u ,9E9M 1 \- (h c \x k \KL)(KL\(s c - sJW) + (h c \(s c - s e ) t \KL) {KL\x k W) 

M-EXMX = i -z — > tj ? (7T \E k B t 0), (51) 

6 2m c ^ M^,/ - Ekl 

where s c and s g are spins of c and c, respectively. The phenomenological approach to the H factor used above 
does not work in the present case since there is no accurate measurement of E1M1 transition rate available as 
input datum to determine the phenomenological parameter so far. Fortunately, evaluation of this special H factor 
from QCD turns out to be easy. Since tt° is a pseudoscalar, the H factor (7r°|i?fe-Bi|0) is nonvanishing only when 
E k E>\ = 5 k iE ■ B/3, and E ■ B is related to the axial- vector anomaly. Therefore 

(n°(ri)\a a E%Bt\Q ) = ^(^WK^ ■ B a |0) = ^(^(^F^F^lO}, (52) 
and the last matrix element can be evaluated by using the Gross- Treiman-Wilczeck formula [34] which leads to 

(AasF^nO) ^ ^^^Uml ( V \a s F^\0) = ^Um% (53) 
V 2 m d + m u ^ V6 

in which the factor (m,d — m u )/(md + m u ) reflects the violation of isospin. To predict the transition rate with these 

g 2 g 2 

expressions, we should determine the relations between the effective coupling constants = — , an = — — and 

47T 47T 

the coupling constant a s appearing in Eqs. (52) and (53). With certain approximations, we can calculate the 
transition rates T{ip' — > J/ip irir) and T(ip' — » J/ip rj) expressed in terms of ole and oim [7, 14], so that cue and 
<xm can be determined by taking the input data (24) and (27) [61]. The determined <xe is approximately 

a E w 0.6, (54) 

while the determination of ctu is quite uncertain because the approximation used in calculating T(\p' — » J /ip rj) 
is rather crude [14]. So we take a possible range [14] 

«e < ctM < 3ae (55) 

to estimate the rate. Since the value of cie in (54) is just about the commonly estimated value of the strong 
coupling constant a s at the light hadron scale, we simply take a s w a E - In this spirit, taking account of the S-D 
mixing (40) in ip', the transition rate of (50) is 

2 



]?(</>' - h c ir°) 



143m^ V aE 



x ■ 



m d i'"u p 

; U m 

m d + m u 



cos6( + Chj - a/2 sin ^./^ + ffi n 

q v \. (56) 



Here we have neglected the state mixing effect in h c which is small [26] since h c is not close to the DD threshold. 
Numerical result in the CK potential model is [24] 

r(V>' ^h c Tr a ) = 0.06(—) keV, BN;' -» h c 7r°) = (2.2 ± 0.2) x 10~ 4 . (57) 

The calculation shows that the dependence of the transition rate on the potential model is mild. 

We know that ir° decays 99% into two photons. Thus the signal in (50) is tp' — ► h c -f-f with M 77 = m^o . If the 
momenta of the two photons can be measured with sufficient accuracy, one can look for the monotonic M 11 as 
the signal. From the branching ratio in (57), we see that, taking account of a 10% detection efficiency, hundreds 
of signal events can be observed for an accumulation of 10 millions of tp' . The backgrounds are shown to be 
cither small or can be clearly excluded [24]. Once the two photon energies Wi and u 2 are measured, the h c mass 
can be extracted from the relation M% — M^, + m^. — 2M^i{u>\ + oj 2 )- 

To have a clearer signal, one can further look at the decay product of h c . It has been shown that the main 
decay channel of h c is h c — ► r] c -f [24]. So that the easiest signal is tp' — > h c n° — > fy c 777- The branching ratio 
B(h c — * ?7 C 7) depends on the hadronic width of h c . In Ref. [24], the hadronic width of h c was studied both in 
the conventional perturbative QCD (PQCD) and in nonrclativistic QCD (NRQCD) approaches. 

We first look at the PQCD result. With the hadronic width obtained from PQCD, Ref. [24] predicts 

B(h c ^n cl ) = (88 ±2)% 

Combining (57) and (58) with the possible range (55) of the undetermined parameter cxm/me, we obtain 

PQCD: B(tp' — > h c ir°) x B(h c — > n c j) = (1.9 — 5.8) x 10~ 4 . (59) 

Signals considering the exclusive hadronic decays modes of rj c are also studied in Ref. [24] . 

Recently, CLEO-c has found the h c state via the channel tp' — > h c n° — * ?7c777 [16, 35]. The measured resonance 
mass is M hc = 3524.4 ± 0.6 ± 0.4 MeV [16, 35] which is consistent with the E835 result at the la level. The 
measured B{ip' — > /i c 7r°) x B(h c ^> r] c j) is [16, 35] 

CLEO-c: B(tp' -» ^ c tt°) x B(^ c -» 77 c7 ) = (3.5 ± 1.0 ±0.7) x 10" 4 (60) 

which is in good agreement with the above theoretically predicted range (59). Future improved measurement 
with higher precision can serve as an input to determine the unknown parameter olmIole- 

NRQCD predicts a larger hadronic width of h c , so it predicts a smaller branching ratio of h c — ► rj c ^j, say 
B(h c -» r/ c7 ) = (41 ± 3)% [24] which leads to 

NRQCD : B(ip' -» /i c tt°) x B(ft c -» r/ c7 ) = (0.9 - 2.7) x 10~ 4 . (61) 

This is also consistent with the CLEO-c result (60) to the present precision. We expect future CLEO-c experiment 
with higher precision to test the PQCD and NRQCD approaches. Since the H factor (53) in this tp' — > h c n° 
process is obtained from the Gross- Treiman-Wilczek relation without taking approximations, the agreement 
between (59) and (60) implies that the above theoretical approach to the MGE factor in (51) is quite reasonable. 
CLEO-c has also studied some exclusive hadronic channels [16, 35]. More accurate measurement of the branching 
ratios of these exclusive hadronic channels may also be compared with the corresponding predictions in Ref. [24] 
to test PQCD and NRQCD approaches. 

IV. NONRELATIVISTIC COUPLED-CHANNEL APPROACH TO HADRONIC TRANSITIONS 

We know that an excited heavy quarkonium state lying above the open heavy flavor threshold can decay into 
a pair of heavy flavor mesons V and V (V stands for the D mesons if the heavy quark is c, and stands for the 
B mesons if the heavy quark is b). This means that there must exist couplings between <£>, V and V shown in 
FIG. 2. With such couplings taken into account, a complete theory of heavy quarkonia satisfying the requirement 
of unitarity should include not only the theory describing the discrete states $, but also the theory describing 
the continuous sector VV as well. Such a theory is the so-called coupled- channel theory. 



(58) 



FIG. 2: Coupling of the heavy quarkonium <E> to its decay channel DT). 

It is hard to study the &-T>-T> vertex shown in FIG. 2 from the first principles of QCD since it is an interaction 
vertex between three bound states. There are various models describing coupled-channel effects, and the two well 
accepted models are the Cornell coupled-channel model [12, 36] and the UQM [26] mentioned in Sec. III. The 
§-D-T> vertex in the UQM is taken to be the 3 Po quark-pair-creation (QPC) mechanism [37], i.e., the creation 
of light quark pair qq is supposed to have the vacuum quantum numbers J PC — ++ ( 3 Pq), and the vertex in 
FIG. 2 is described by the 3 Po sector of the overlapping integral between the three bound-state wave functions 
with an almost universal coupling constant 'Jqpc ~ 3.03 [37]. The parameters in the UQM are carefully adjusted 
so that the model gives good fit to the cc and bb spectra, leptonic widths, etc. It has been shown that the QPC 
model even also gives not bad results for OZI-allowed productions of light mesons [37, 38], which will be relevant 
in the calculation of the hadronic transition amplitudes in FIG. 3(e) and 3(f) in the coupled-channel theory. So 
we take the UQM in this section. 
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FIG. 3: Diagrams for hadronic transitions in the coupled- channel approach, (e) and (f) are new transitions mechanisms beyond the 
QCD multipole expansion. Quoted from Ref. [39]. 



In the UQM, the whole Hilbert space is divided into two sectors: namely, the confined sector |3>o! A) labelled 
by the discrete quantum number A (say a, n, L, J) and the continuous sector \T>T>; v) labelled by the continuous 
quantum number v (say the momentum). The state |<!>o; A) is just the eigenstate of the Hamiltonian Hq in the 
naive single-channel theory with the eigenvalue (the bare mass); i.e., 

H |$ ;A) =M A ° |$ ;A), (62) 

and the state \T>T>; v) is a state with two freely moving mesons T> and T>, which is the eigenstate of the kinetic- 
energy Hamiltonian Hq with the energy eigenvalue E v ; i.e., 

H c q\VV;v) =E v \VV-v). (63) 

The total Hamiltonian H of the system contains Hq, Hq and the quark-pair-creation Hamiltonian Hqpc which 
determines the OZI-allowed $-2?-2? vertex and mixes the two sectors. H can be written as 

H 0\ / Ht 



H =Wh S ) + W qpc T )■ <«> 

where the first and second rows stand for the confined channel and the continuous channel, respectively. Note 
that 

($05 MHqpc\*o; A') = 0, (VV- y\H QPC \VV- v 1 ) = 0. (65) 




FIG. 4: The self-energy Tlx\' from the T> meson loop. 

With Hqpc introduced, there will be a self-energy n^>,' of the quarkonium <I>o contributed by virtual loops 
of T> mesons. This is shown in FIG. 4. The self-energy is not necessarily diagonal, i.e., A and A' may be 
different. This causes the state mixings. For states below the threshold, the self-energy n^y is 

n f {^\\H Q p C \VV; V )(VV- lV \H Q pc\^\') , 

HAA ' =-J M^E V ^ (66) 

Now the total mass matrix of the quarkonium state is 

Af A v = ^AV+n A v. (67) 

Let ct\\i be the matrix diagonalizing M\\i, and M\ be the diagonal matrix element. The physical quarkonium 
state |$; A) is the eigenstate of H with the energy eigenvalue My, i.e., 

F|$;A)=M A |$;A). (68) 

The eigenstate |$; A) can be expressed as a superposition of |<!>o; A) and \T>V; v): 

|$;A) =^a A v|*o;A') + / c x ( V )\VV; V ) dv, (69) 



in which all possible open heavy flavor mesons T> (T>) composed of the heavy quark Q and all possible light 
quarks q (q) should be included. 



The state mixing coefficient aw is related to a xx * by [26] 



a\\> = N x a xx ,, 



(70) 



where [26] 



1 + 



/ 



{Vf)-v\H Q pc\§ ;\) 



M x -E v 



1 -1/2 



dv 



(71) 



is a normalization coefficient, and iV| determines the probability of finding the confined sector |$ ;-^) m the 
physical state |$; A). The calculation of a's and N\$ is tedious, and the results for various cc and bb states are 
given in Ref. [26]. The mixing coefficient c\(v) is related to a\\> by [26] 



caM = Y i a xx ,(VV;v\H QPC \$ ;\')/(M x - E v ). 



x 



(72) 



In the single-channel approach, the energy eigenvalues M° of high lying quarkonium states predicted by 
potential models are usually higher than the experimental values. In the coupled-channel theory, the self-energy 
Ha a' usually causes M\ < M°. Since coupled-channel corrections M\ — M" arc unimportant for states lying 
much lower than the VV threshold but are relatively important for states close and above the VV threshold, 
coupled-channel theory does improve the prediction for the energy spectra. The UQM coupled-channel theory 
has been applied to obtain successful results of heavy quarkonium spectra, leptonic widths, etc., for the cc and 
bb systems [26]. 

The formulation of the theory of hadronic transitions in the framework of the UQM was given in Ref. [39] . 
Let H pair = Hqpc + Hq PC , H = H Q for the confining sector and H = H^ for the continuous sector. In the 
framework of UQM, the S matrix element (15) becomes [39] 



(F\S\I) 



-i2nS(E F 
x (H 2 + H pair 



1 



(H 2 + H pa i r ) 



1 



E T -H + td - Hi 



(H 2 + Hpair) 



E!-H + id - Hi 



(if 2 + H palr ) 



(73) 



For isospin-conserving 7T7t transitions (dominated by E1E1 gluon emissions), we take the electric dipole term in 
H 2 [cf. Eq. (11)]. Note that, in Eq. (73), the creation of the two pions can come not only from the conversion of 
the two emitted gluons (OZI-forbidden mechanism) via the two H 2 , s, but also from the OZI-allowed mechanism 
(VDn; v\Hqpc\DT>\ v) directly from the light quark lines [cf. FIG. 3(e)-3(f)]. Note that the two gluons can only 
convert into two pions (not one pion) due to isospin conservation. Thus these two pion creation mechanisms 
contribute separately. The mr transition S matrix element between two physical quarkonium states |$;A/) and 
|$; X F ) is [39] 



($;A F ;7r(fc 1 )7r(fc 2 )|5|$;A / ) = -i2nS{M f + E ni + E V2 - Mj) ^ a Xl y a XF \> f ( $o; A^; 7r(fei)7r(fe 2 ) 
1 „ „+ 1 „ 1 



Ho 



'Mj-Ho + ido-Hi 



-H 2 + H} 



QPC 



-Ho 



-Ho 



1 



+H 2 - 



1 



-H. 



t 

QPC 



Mi -H + id - Hi Mi -H + id - Hi Mi - H 
Hi — — Hq pc 



Hqpc 



Mt - Ho + id - Hi Mi — H + ido — Hi Mi - H 
+H Qpc— . ~ ~ H2 ~ r \ . ~ —Hqpc 



+H, 



QPC 



Mi -H a + ido -Hi Mi-H + id a - Hi 

Hqpc— 77 Hqpc— jtHqpc 



Mj-Hq 



Mi -Ho 



Mi -Ho 



M!-H + id - Hi 



-H 2 



(74) 



where E 7T1 and E 7T2 are energies of the two pions. The Feynman diagrams corresponding to the terms in (74) 
are shown in FIG. 3 in whic FIG. 3(a)— Fig. 3(d) are diagrams corresponding to the first four terms in (74), and 
FIG. 3(e)— FIG. 3(f) are diagrams for the last term in (74). For convenience, we shall call the first four terms in 
(74) the MGE part, and call the last term in (74) the quark-pair-creation (QPC) part. 



We see that (74) contains much more channels of tttt transitions than the single-channel theory does. In the 
MGE part, FIG. 3(a) is similar to FIG. 1 but with state mixings, so that the single-channel amplitude mentioned 
in Sec. IIIA is only a part of the first term in (74). In the QPC part, the last term in (74) is a new pion creation 
mechanism through Hqpc irrelevant to MGE. Thus in the coupled-channel theory, tttt transitions between heavy 
quarkonium states are not merely described by QCDME. 

Since the state mixings and the QPC vertices depending on the bound-state wave functions are all different 
in the cc and the bb systems, the predictions for T(T' — > T tttt), T(T" — > T tttt) and r(T" — > T' tttt) by taking 
T expt {tjy — > J/ip tttt) as input will be different from those in the single-channel theory. Such predictions were 
studied in Ref. [39] in which the same potential model as in Ref. [26] is taken for avoiding the tedious calculation 
of a's and N\s. Note that for a given QPC model, the QPC part in (74) is fixed, while the MGE part still 
contains an unknown parameter C\ in its hadronization factor after taking the approximation (36). Since there 
is interference between the MGE part and the QPC part, the phase of C\ will affect the result. Let 

C 1 =\C 1 \e i ». (75) 

Two input data are thus needed to determine |Ci| and §. In Ref. [39], the data of the transition rate and 
M„„ distribution in tp' — > J/tp tttt are taken as the inputs. Considering the experimental errors in the M„ n 
distribution, ■& is restricted in the range — 1 < cos?? < —0.676. The details of the calculation are given in Ref. 
[39] in which the V meson states D (B), D* (B*), and D** (£?**) are taken into account. The so predicted 
transition rates T(T' — * T tttt), T(T" — > T tttt), and r(T" — * T' tttt) for cosi9 = — 1 and cos'd — —0.676 are 
listed in TABLE IV together with the updated experimental results for comparison. 

TABLE IV: The predicted rates T(T' -> X tttt), r(T" -> X tttt), and T(X" -> X' tttt) ( in kcV) in the coupled-channel theory with 
cosi? = —1 and cos$ = —0.676. The corresponding updated experimental values of the transition rates quoted from Ref. [11] are 
also listed for comparison. 
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1.26 ± 0.40 



We see that the obtained T(T' — > T tttt) is in good agreement with the experiment, and the results of L(T 
T tttt) and r(T" — > T' tttt) are in agreement with the experiments at the level of 2a and 2.4c, respectively. 
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FIG. 5: Comparison of the coupled-channel theory predicted curve of dV(T' — > X tttt) / 'dM-^n with the ARGUS data [40]. The solid 
and dashed-dotted lines stand for cos 9 = —1 and cos6> = —0.676, respectively. The dashed line is the naive single-channel result for 
comparison. Quoted from Ref. [39]. 

Next we look at the predicted M TI distributions. It is pointed out in Ref. [40] that there is a tiny difference 
between the the measured M n7r distributions in iji 1 — ► J ftp tttt and T' — > T tttt. In the single-channel theory, the 
formulas for these M M distributions are the same with the same value of |Ci|. Ref. [40] tried to explain the 
tiny difference by taking the approach to the H factor given in Ref. [4] in which there is a parameter k which 



is supposed to run. However, the running of k is not known theoretically, so that it is not clear whether the 
running of K from the scale Afy — Mj/^ — 590 MeV to the scale M-p — My = 560 MeV can really explain the 
tiny difference or not. Furthermore, as we have seen in Sec. IIIC that the approach given in Ref. [4] is ruled out 
by the recent BES and CLEO-c experiments. In the present coupled-channel theory, once the values of \C\\ and 
■& are determined by the input data of ip' — > J/ip tttt, the M„ v distribution of Y' — > T tttt is definitely predicted. 
The comparison of the predicted distribution with the experimental data given in Ref. [40] is shown in FIG. 5. 
We see that the agreement is good, so that the coupled-channel theory successfully predicts the tiny difference. 

However, the situation of the M T7r distributions of T" — > T tt + it~ and T" — > T' tt + it~ are more complicated. 
The single- channel theory predicts distributions similar to FIG. 5 for these two process, i.e., the distributions 
are peaked at the large M W7r region. The CLEO data shows a clear double-peaked shape for the M nit distribution 
of T" — > T 7r + 7r~ [cf. Fig. 6(a)] [41, 42]. The coupled-channel theory does enhance the low-M T7r region a little, 
but is far from giving a double-peaked shape as is shown by the solid and dashed-dottd curves in FIG. 6(a). 
Actually, this situation is not only for the coupled-channel theory based on the UQM. The Cornell coupled-channel 
model is not substantially different from the UQM [39] . Compared with the UQM, the Cornell coupled-channel 
model leads to relatively larger S — S mixings but smaller S — D mixings after taking the same experimental 
inputs. So that the Cornell coupled-channel model gives even smaller enhancement in the low-M^ region. Thus 
the transition T" — > T it + ir~ needs further investigations with new ideas although the predicted transition rate 
F(T" — > T 7r + 7r~) is consistent with the CLEO data at the 2a level. 




FIG. 6: Comparison of the coupled-channel theory predicted curves of gT(T" — * T tt + tt )/dM 1TW [FIG. 6(a)] and dr(T" — > 
T' 7r+7r-)/dA/ w7r [FIG. 6(b)] with the CLEO data [41]. The solid and dashed-dotted lines stand for costf = -1 and cos0 = -0.676, 
respectively. The dashed line is the naive single-channel result for comparison. Quoted from Ref. [39]. 

There have been various attempts to explain the double-peaked shape. Ref. [43] assumed the existence of a 
four-quark state Ti having nearly the same mass as T" and coupling strongly to T"ir and Ttt, and the dominant 
transition mechanism is suggested to be T" — > Ti + tt —> T + tt + tt which enhances the low-M W7r distribution. The 
branching ratio of T(4S f ) — > Ti + tt is estimated to be roughly 1%, so that the assumption can be experimentally 
tested by searching for the Ti state in T(4S') decays. This idea was carefully studied in Ref. [44] taking account 
of the final-state tttt interactions and got a double-peaked shape but the low--M 7r7r peak is not at the desired 
position. A slightly modified model of this kind was proposed in Ref. [45]. So far the assumed four-quark 
state is not found experimentally. Another attempt was made in Ref. [46] assuming that the coupled-channel 
contributions are strong enough in T" — ► T tttt that there is a considerably large QPC part in the transition 
amplitude, and the interference between it and the MGE part may form a double-peaked shape by adjusting 
the strength of the QPC part. However, as we mentioned above that the strength of the QPC part is fixed once 
a QPC model is given, and the systematic calculation in Ref. [39] shows that the QPC part is actually much 
smaller than what was expected in Ref. [46]. Recently, attempts to explain the double-peaked shape by certain 
models for a light a meson resonance at around 500 MeV in the final state tttt interactions with [47] and without 
[48] using the Breit-Wigner formula have been proposed. By adjusting the free parameters in the models, the 
CLEO data on the M„ distributions in T" — > T tttt and T" — > T' tttt can be fitted. However, the models need to 



be tested in other processes. Therefore, the T" — > T 7T7t transition is still an interesting process needing further 
investigations. 

We would like to mention that the calculations mentioned above concerns the wave functions of some excited 
states of heavy quarkonia, the heavy flavored mesons T>, and the pions. Nonrelativistic potential model cal- 
culations of these wave functions may not be so good. Therefore the nonrelativistic coupled-channel theory of 
hadronic transitions in Ref. [39] still needs further improvements. 



V. APPLICATION OF QCD MULTIPOLE EXPANSION TO RADIATIVE DECAYS OF J/ip 

In the preceding sections, QCD multipole expansion is applied to various hadronic transition processes in which 
the initial- and final-state quarkonia $/ and &p are composed of the same heavy quarks. In this case, the dressed 
(constituent) quark field *£(x,i) does not actually need to be quantized. Now we generalize the QCDME theory 
to processes including heavy quark flavor changing and heavy quark pair annihilation or creation. Then the 
quantization of the ^f(x.t) is needed. This has been studied in Ref. [6], and the obtained canonical commutation 
relation is [6] 



[V(x,t),&(x',t)] =5 3 (x-x'). 
To include the electromagnetic and weak interactions, we generalize the Hamiltonian as 



II — Hqq D -Hint j 
r(l) 



Hint — H ' QCD + H em + Htfr, 



(76) 

(77) 
(78) 



in which Hq^ d and Hq CD are defined in Eqs. (9) and (11), and 



r(i) 



Hp 



Hw — 



e / d 3 x ^(x,t)^QAJx,t)^{x,t), 
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t 3 -sm z 6 w Q)Z ll {x)^{x 1 t) 



(79) 



where Q is the electric charge operator of the heavy quark, A u is the photon field, g and ti are, respectively, the 
weak SU(2) coupling constant and generator, 9w is the Weinberg angle, and e = gsin9\y is the electromagnetic 
coupling constant. 




J/i, 



(a) 



J/t/> 




(10 
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FIG. 7: Feynman diagrams for the radiative decay process J/ip 
are all bound states. 



■7 + 17. The intermediate states between two vertical dotted lines 



Let us take the application of the generalized theory to the radiative decay process J /if) — * 7 + 77 as an example. 
This process has been studied in the framework of perturbative QCD and nonrelativistic quark model in Ref. 
[49], but the predicted rate is significantly smaller than the experimental value. We know that the momentum 
of the 77 meson in this process is q v = {M^,^ — m^)/(2Mj/^) = 1.5 GeV. Suppose the r\ meson is converted from 
two emitted gluons from the heavy quark. The typical momentum of a gluon is then k ~ q-q/2 ~ 750 MeV. This 
is the momentum scale that perturbative QCD does not work well but QCDME works [6]. So we can calculated 



the rate of this decay process using QCDME. The Feynman diagrams for this process are shown in FIG. 7, in 
which the intermediate states marked between two vertical dotted lines are all treated as bound states in this 
approach. In this sense this approach is nonperturbative, and also for this reason the contributions of the three 
diagrams in FIG. 7 are different. 

In QCDME, this process is dominated by the E1M2 gluon emissions. So the H factor (conversion of the two 
gluons into 77) is 

gEQM^D.B^), (80) 
where Dj = dj — g s (X a /2)Aj is the covariant derivative. The operator in (80) can be written as 

EJDjB? = dj(E^Bt) (DjEj)Bf. 
It is argued in Rcf. [4] that the second term is smaller than the first term, so they suggested the approximation 

g E g M (v\E«D 3 Bt\Q) c iq^E^B^O). (81) 
This matrix element can then be evaluated by using the Gross- Treiman-Wilczeck formula [34] and we obtain [6] 

g E g M (v\E^D,Bt\0) ~ i 9 E ^ M ^- g^Am^cosgp - \/2sin0p)*y, (82) 
where Op is the mixing angle in the pseudoscalar nonet, i.e., 

77 = i]s cosdp — 771 sin#p, 7/ = i]$ sin Op + 771 cos Op, (83) 

with 

1 1 
m = —i=(ui5U + drf 5 d - 2s7 5 s), 771 = —=(wy 5 u + d-f 5 d + sj 5 s). (84) 
V6 V3 

As in Sec. HID, we take a s ~ a E = 0.6. It is shown in Ref. [6] that the contribution of FIG. 7(a) is larger 
than those of FIG. 7(b) and FIG. 7(c). Thus as an approximation, we only consider the main contribution of 
FIG. 7(a). 



The calculated J ftp — > 777 rate is [6] 

1 (a M \ \q n ? ( 2eQ 



r(J/r/> -> JV) y aE j y 3 ^/g TOc 



4tt 2 n2 
-^-y=Um*(cos P - V2sin6»p)^/i^ 



(85) 



where 



h LPjP ST ( R nF\ r \ R Kl)( R KL\ r \ R nih) f /px fRR s 

in which Rn^, R n i and R' KL are radial wave functions of the initial-, final-, and intermediate-quarkonium states 
in FIG. 7(a), respectively. / n ;(0) is the wave function at the origin of the final-quarkonium state. We shall take 
into account the first five terms in the summation ^""^ ^ionO- As ^ s well-known that /„o(0) can be determined 

n 

by the datum of the related leptonic width T(ijj(n 3 Si) — ► e + e~). For n — 1,2, the so determined /io(0) and 
/2o(0) are smaller than the ones predicted by the Cornell potential model by almost the same factor 0.57 [12]. 
It is expected that this discrepancy may be explained by QCD corrections. For n > 3, the states are above the 
threshold and state mixings will be significant, so that the data of T(ip(n 3 Si) — > e+e~) are not useful. We expect 
that QCD corrections will not vary seriously with n as is inspired by the cases of n = 1,2. Then we can calculate 
/io(0),- • -/5o(0) using the Cornell potential model, and then multiply the obtained results by the same factor 
0.57 to obtain the correct values of them. 

The factor (cos Op — v^sinflp) in (82) concerns the effective r\-g-g vertex in the hadronization gg — > 77. This 
is somewhat similar to the effective 77-7-7 vertex in 77 — > 77. We may take the determined value of Op from the 
77 — > 77 and 77' — + 77 data, which is Op w —20° [11]. With this value of Op, we get 



T(J/i> -» 777) = 0.041 0^) keV - ( 87 ) 



With the experimental datum T to t{J/ip) — 91.0 ± 3.2 keV [11], we obtain 



B( J/ip -» 777) = (4.5 ± 0.2) x 10" 



4/ OLM 
OLE 



(88) 



The experimental value of this branching ratio is [11] 



(8.6 ±0.8) x 10" 



expt 



We see that for olm/ole ~ 1.9, the predicted branching ratio agrees with the experimental value. 

Note that the value of um/q-e and 9p are not so certain, and we do not know how good the approximation 
(81) really is. To avoid these uncertainties, we can take the ratio of T{J/ip — > 777) to another E1M2 transition 
rate T(ip' — > J/iprj). The theoretical prediction is [6] 



_ r(J/^ 7 r?) 

" - r(v - J/Vr?) 



-(eQ)^(JA^ 7^)17^ 



Efclll 
"lOnO 



0.012. 



(89) 



/ 2010 l 



In this ratio, the uncertainties mentioned above are all cancelled, so that R v just tests the MGE mechanism in 
this approach. The corresponding experimental value is [11] 



R, 



0.009 ±0.003. 



(90) 



expt 



We see that the agreement is at the la level. Since we have seen in Sec. HID that the calculation of the MGE 
factor mentioned in Sec. Ill is quite reasonable, the agreement of (89) with (90) implies that MGE mechanism 
for this radiative decay process is also reasonable. 

The above approach can also be applied to the radiative decay process J/ip — > 77/ . From (83) we see that the 
J/ip — > 77/ decay rate is 



vim - W) - (m] ^ (-^- 



4tt 2 
3V6 



Urn*, {V2 cos Bp + sin 6 P ) ^ h 1 ^ 



(91) 



Since there is no ip' — ► J/iprj' available (not enough phase space), we cannot have a ratio similar to R n which 
exactly tests the MGE mechanism. We can define the ratio 



q(J/V> -> 777') 



q(J/V> -> 77?)) 



7 2 ,(\/2cosflp + sin 6» P ) 



i 2 (cos6*p — \/2sin0p) 



(92) 



Taking 9 P w —20° determined from the 77 — > 77 and 77' — > 77 rates, we predict 



i? v = 0-044. 



The corresponding experimental value of R v i is [11] 



0.044 ±0.010. 



expt 



(93) 



(94) 



We see that this prediction is also in agreement with the experiment. 

It has been shown in Rcf. [6] that the contribution of the above mechanism to the isospin violating radiative 
decay J/ip — > 7 + ir° is negligibly small. There is another important mechanism giving the main contribution to 
J/ip —* 7 + 7T°. It is the p° meson dominance mechanism J/ip — > p * + 7r° — » 7 + 7r°. This has been studied in 
Ref. [58] , and the result is close to the experimental value [62] . 

We would like to mention that this approach is not suitable for T — > 777 since the typical gluon momentum in 
this process is k ~ q v /2 <~ 2.4 GeV at which perturbative QCD works, while QCD multipole expansion does not. 
Studies of the processes T — > 777, T — > 777' and T — ► 7/2 (1270) have been carried out in Ref. [50]. Application 



of this approach to ip' — > 777 is more complicated since both relativistic and coupled-channel corrections are 
important in this process. Thus developing a relativistic coupled-channel is desired. 

QCDME can also be applied to study the direct photon spectrum in J/ip — ► 7 + hadrons near x ~ 1, where 
x = 2u 7 /M with w 7 the energy of the photon and M the mass of the quarkonium. Conventional study of the 
process J/ip(T) — > 7 + hadrons are based on perturbative QCD calculation of J/ip{T) — > 7 + 5 + g in the Born 
approximation [51]. The direct photon spectrum is expressed as 

1 dT(J/ip(T) -> 7 + hadrons) 
r to t cfe 

For T, the obtained result is in good agreement with the experiment [52-54], while for J/ip the obtained dis- 
tribution is too hard, i.e., in the range x > 0.8, the obtained distribution is much larger than the experimental 
values [55]. In J/ip — > 7 + hadrons, the typical gluon momentum at x ~ 1 is fc ~ 770 MeV, so that we can apply 
QCDME to it for x > 0.9. The calculation was done in Ref. [56], and the obtained direct photon spectrum for 
x > 0.9 is very close to the experimental values [56]. For x < 0.9, the typical gluon momentum is too large for 
QCDME to work. A successful theory for the whole range of x is still expected. 



VI. SUMMARY AND OUTLOOK 



In this paper, we have reviewed the theory and applications of QCDME. We see from Sees. Ill— V that 
nonrelativistic QCDME theory gives many successful predictions for hadronic transition and some radiative 
decay rates in heavy quarkonium systems. Even the simple nonrelativistic single-channel theory can work well 
for many processes. Although the single-channel approach gives too small rates for T' — > Ttttt, T" — > Ttttt, and 
T" — ► T'tttt (cf. TABLE I), nonrelativistic coupled-channel theory improves the prediction (cf. TABLE IV). 
We summarize the above mentioned main successful predictions for the transition and decay rates in TABLE V 
together with the corresponding experimental results for comparison. 

TABLE V: Summary of the predictions for transition and decay rates in the nonrelativistic QCD multipole expansion 
approach together with the corresponding experimental results for comparison. 



Theoretical predictions 



Experimental data 



Places in the text 



R' 



R" 



T(T' -> T 7T7T) 
T(T" T Tin) 
T(T" -> T' tttt) 
T(T' ->• T 77) 

r(T" -> t n) 

\f^o(bb)/™ b \ 2 \q(bb)\ 3 

l/2 1 01o( C c)/™cf 

\ti£o(bb)/m b \ 



|g(cg)|3 
\q(bb)\ 3 



r( X 6(2 3 P ) - X6(1 3 P ) 
r( X 6(2 3 P ) - X6(1 3 P 2 ) 777T 
r( X 6(2 3 P!) - X b(l 3 Pl) ™ 
r( X 6(2 3 Pl) -> X6(l 3 -P 2 ) 7V7V 
r( X 6(2 3 P 2 ) -» X6(1 3 P 2 ) 7V7V 

r(t/>(3770) -» J/ip 



B(ip' h c -K°)B{h c -» r?c7) 
R v = Y{J^ - 7 r?)/r(V - J/V r?) 

j? v = r(j^ - 7 r?')/r(V' - v) 



13 keV 

1.0 keV 

0.3 keV 
0.022 keV 
0.011 keV 

0.0025 

0.0013 

0.4 keV 
0.002 - 0.02 keV 

0.4 keV 
0.001 - 0.01 keV 

0.4 keV 
(32 - 147) keV 

(1.9-5.8) x 10~ 4 
0.012 
0.044 



12.0 ± 1.8 keV (PDG) 
1.72 ±0.35 keV (PDG) 
1.26 ±0.40 keV (PDG) 

< 0.086 keV (PDG) 

< 0.058 keV (PDG) 

< 0.0098 (BES, PDG) 

< 0.0065 (BES, PDG) 



(0.83 ± 0.22 ± 0.08 ± 0.19) keV 

(0.83 ± 0.22 ± 0.08 ± 0.19) keV 
(80 ± 32 ± 21) keV (BES) 
(50.5 ± 16.9) keV (CLEO-c) 
(4 ±0.8 ±0.7) x HT 4 (CLEO-c) 
0.009 ± 0.003 
0.044 ± 0.010 



TABLE IV 
TABLE IV 
TABLE IV 
Eqs.(29)(30) 
Eqs.(29)(30) 

Eqs.(32)(31) 

Eqs.(32)(31) 

TABLE II 
TABLE II 
TABLE II 
TABLE II 
TABLE II 
TABLE III,Eq.(45) 
Eq.(48) 
Eqs.(59)(60) 
Eqs.(89)(90) 
Eqs.(92)(94) 



In addition, the prediction for the M vn distribution in T' — > T7T7T from the nonrelativistic coupled-channel 
theory is in good agreement with the data (cf. FIG. 5). 

However, despite of the above success, there are experimental results which this simple nonrelativistic approach 
cannot explain. CLEO experiment shows a clear double-peak shape for the M n7r distribution in T" — > T 7T7t [cf. 



FIG. 6(a)]. Nonrelativistic coupled-channel correction is so small that it cannot account for this shape. Whether 
this double-peak shape can be explained by the final state tttt interactions or it is caused by other physical effects 
is still not clear yet. Further investigation is needed. 

Another problem is that, in the nonrelativistic single-channel approach, the S- wave to S- wave transitions rates 
are contributed only by the C\ term in Eq. (19) due to the orbital angular momentum selection rule, i.e., the 
obtained tttt angular correlation is isotropic in the laboratory frame. However, experiments on — ► J/tp tttt 
[57] shows a small angular dependence (a 0.2% ingredient of D-wave) of the tttt angular correlation which cannot 
be explained by the nonrelativistic single-channel theory [63]. Theoretically, the angular dependence of the 
transition rates may come from: (a) coupled-channel corrections [state-mixing leads to the C2 term (D-wave) 
contributions], and (b) relativistic corrections [orbital angular momentum no longer conserves in the relativistic 
theory]. Actually, the sizes of the corrections (a) and (b) are of the same order of magnitude, and there is 
interference between them. Therefore to obtain a theoretical prediction for the 7T7t angular correlation, both 
(a) and (b) corrections should be taken into account. This means that a systematic relativistic coupled-channel 
theory of hadronic transitions is expected. So far there is still no such a theory due to the difficulty of dealing 
with the two-body bound-state equation in relativistic quantum mechanics. There have been various attempts 
to solve the relativistic two-body problem. Making effort on developing a systematic relativistic coupled-channel 
theory of hadronic transition is really important. 

In a word the nonrelativistic theory of QCDME approach is not the end of the story. Further development is 
needed. 
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